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ABSTRACT 

We explore the mid-infrared (mid-IR) through ultraviolet (UV) spectral energy distributions (SEDs) 
of 119,652 luminous broad-lined quasars with 0.064 < z < 5.46 using mid-IR data from Spitzer 
and WISE, near-infrared data from Two Micron All Sky Survey and UKIDSS, optical data from 
Sloan Digital Sky Survey, and UV data from Galaxy Evolution Explorer. The mean SED requires 
a bolometric correction (relative to 2500 A) of BC '— 2.75 ± 0.40 using the integrated light from 

1 /im-2 keV, and we further explore the range of bolometric corrections exhibited by individual objects. 
In addition, we investigate the dependence of the mean SED on various parameters, particularly the 
UV luminosity for quasars with 0.5 ^ ^ ^ 3 and the properties of the UV emission lines for quasars 
with z > 1.6; the latter is a possible indicator of the strength of the accretion disk wind, which is 
expected to be SED dependent. Luminosity-dependent mean SEDs show that, relative to the high- 
luminosity SED, low- luminosity SEDs exhibit a harder (bluer) far-UV spectral slope (auv)> a redder 
optical continuum, and less hot dust. Mean SEDs constructed instead as a function of UV emission line 
properties reveal changes that are consistent with known Principal Component Analysis (PGA) trends. 
A potentially important contribution to the bolometric correction is the unseen extream-UV (EUV) 
continuum. Our work suggests that lower-luminosity quasars and/or quasars with disk-dominated 
broad emission lines may require an extra continuum component in the EUV that is not present (or 
much weaker) in high-luminosity quasars with strong accretion disk winds. As such, we consider 
four possible models and explore the resulting bolometric corrections. Understanding these various 
SED-dcpendcnt effects will be important for accurate determination of quasar accretion rates. 
Subject headings: catalogs — infrared: galaxies — methods: statistical — quasars: general 
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1. INTRODUCTION 

Most massive (bulge-dominated) galaxies are believed 
to harbor a supermassive bl ack hole at their centers 
(jKormendy fc RichstonellT995l ). For the most part, these 
black holes are passive, but when the host galaxy's 
gas loses angular momentum it acc retes onto the blac k 
hole, creating a luminous quasar (|Lvnden-Belll I1969D . 
This gas can be disrupted i n a number of ways, 
inclu d ing galaxy me r gers (e.g. IKauffmann fc Haehneltl 
20001: iHopkins et al.l l2006f) and secular fueling (e.g. 
Hopkins fc Hernguistl [2009, and references therein). 
Most of the energy output is due to an accretion disk 
that forms about the black hole which either directly 
or indirectly produces radiation across nearly the entire 
electromagnetic spectrum (EM) ([Elvis et al.lll994t ). As a 
result of the quasar process, some of the energy is "fed 
back" into the galaxy, possibly disrupting the galaxy's 
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gas supply (|Silk fc Reeslll998HFabianlll999( ) and quench- 
ing accretion of matter onto the black hole itself (e.g., 
iDi Matteo et al.ll2005HDiorgovski et al.ll2008D . 

Theoretical models of quasar feedback make predic- 
tions that are bas ed on bolometric luminosity (e.g., 
IHopkins et al.|[2006[ ) , while what can be measured is usu- 
ally a monochromatic luminosity, L\. The bolometric 
luminosity, iboi, is given by the integrated area under 
the full spectral energy distribution (SED). The ratio 
between Lboi and L\ defines a "bolometric correction" , 
which can be applied generically if there is reason to 
believe that the quasar SED is well-known (if not well- 
measured for an individual object). Currently the best 
quasar S EDs are bas ed on only te ns or hundreds o f bright 
quasars (jElvis et al. 1994; Richa rds et al.ll2006al ). How- 
ever, over 100,000 luminous broad-lin e quasars have been 
spectroscopically confirmed (Schneider et al.l [20101) . As 
it is dangerous to assume that we can extrapolate the re- 
sults from a few hundred of the brightest quasars to the 
whole population, it is important to expand our knowl- 
edge to cover the ful l quas ar sar nple. 

While lElvis et al.l p99l and lllichards et all ([2006aD 
provide the most complete SEDs in terms of number of 
objects and overall wavelength coverage, further under- 
standing of the SED comes from a variety of other inves- 
tigations from both spectroscopy and multi-wavelength 
imaging across the EM spectrum. For example, compos- 
ite SEDs, from the radio to the X-ray, for 85 optically and 
radio selected bright quasars {\og{vL,y) \ -> 44) 

with z < 1.5 were made by IShang et al.l ([20111 ). and 



the bol ometric corrections for this sample were tabu- 
lated bv lRunnoe et all (|2012[ ). iStern fc Lao3 (|2012[ ) con- 
struct SEDs for over 3500 low-z {z < 0.2) type 1 ac- 
tive galactic nucleus (AGN) with broad Ha lines ranging 
from the near-infrared (near -IR) up to the X- ray with 
log (lyL^) I = > 42, while lAssef et all ()2010f ) present 

A — zoUU A 

SEDs using over 5300 AGNs spanning the mid-infrared 
(mid-IR) up to the X-ray, with redshift going up to ~ 
5.6, from the AGN and Galaxy Evolution Surv ey (AGES; 
iKochaneket al.ll2012D . While in the mid-IR. IPeo et alJ 
(|2011h produced a composite spectrum using 25 luminous 

t ype 1 quasars at z ^ 2. 

iVanden Berk et al.l (pOOlh used 2200 Sloan Digital Sky 
Survey (SDSS) spectra in the redshift range of 0.04 < 
z < 4.79 to construct a mean quasar spectrum cover- 
ing wavelengths from 800A to 8555A in the rest frame. 
From this they found that the mean UV continuum is 
roughly a power-law with a^, — —0.44 (f^ ex v"). To ex- 
plore the fa r -UV ( FUV) region of the quasar spectrum, 
iTelfer et al.l (pOOa ) used - 330 Hubble Space Telescope 
spectra of 184 quasars, with z > 0.33, covering rest 
frame wavelengths from 500A to 1200A, and found an 
anti-correlation between the spectral index of the FUV 
(apuv) and the luminosity at 2500A. This work was later 
supplement e d at fo wer redshift and lower luminosity by 
IScott et all (|2004l ). who used 100 FUSE spectra cover- 
ing the FUV. By combining their data with that from 
iTelfer et al.l (|2002[ ). they found an anti-correlation that 
can be characterized by 



TABLE 1 
Quasar Catalog Format 



apuv = 21.02 -0.49 log 



llOOA llOOA 

erg s~^ 



(1) 



At shorter wa v elengt hs, using 73 quasars, 
lAvni fc Ta nanbaum' (1982) found a dependency of 
the spectral index between the optical and the X- 
ray and a quasar's luminosity. This relation was 
further stud i ed by a number of authors, including 
iSteffen et al.l ()2006[ ). who used 333 quasars with z < 6 
and \og{i'L^) \ °> 42, to find the relationship 

between the UV and X-ray luminosities to be: 



log(L2kev) = (0.721±0.011)log(i^^^^^)+(4.531±0. 

(2) 
while iJust et all ()2007[ ) have extended this res ult to 
higher luminosities. Recently, iLusso et"an (|2010( ) com- 
pleted a similar study using 545 X-ray selected quasars 
and found a similar relationship. In the X-r ay regime, the 
mean SED appears to have ai, ^ —1 fe. g..lGeorge et al" 
20001) before cutting off at - 500 keV (jzdziarski eFaT 



19951 ) 



To improve our understanding of the mean SED (and 
its range), we have constructed quasar SEDs, spanning 
from the mid-IR (~30 /im in the rest frame) to the FUV 
(~300 A in the rest fram e), for quasars cataloged by the 
SDSS (jYorket al.ll2000[) . The data used for this anal- 
ysis are presented in Section [2] Section [3] describes all 
the corrections applied to the data. Section |4] gives an 
overview of our data analysis and construction of mean 
SEDs. Section [5] presents our findings for individual and 
mean bolometric corrections, and our discussion and con- 
clusions are presented in Section [6] and Section [7] respec- 
tively. Throughout this paper we use a ACDM cosmol- 
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Previously published name 

Right ascension in decimal degrees (J2000) 

Declination in decimal degrees (J2000) 
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Redshift 
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BEST SDSS u band PSF magnitude 


6 


Error in u magnitude 
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BEST SDSS g band PSF magnitude 


8 
9 
10 


Error in g magnitude 

BEST SDSS r band PSF magnitude 

Error in r magnitude 



Note. — This tabic is published in its entirety in the 
electronic edition of the online journal. A portion is shown 
here for guidance regarding its form and content. 

ogy with Ho = 71 km s"^ Mpc"\ fl^ = 0.734, and 
Qm = 0.266, consistent with t he Wilkinson Micr owave 
Anisotropy Probe 7 cosmology (jJarosik et al.|[20TTI ). 

2. DATA 

Our sample start s with the SDSS-DR7 quasar catalog 
by [Schneider et all ()2010D . containing 105,783 spectro- 
scopically confirmed broad-lined quasars. Each quasar 
has photometric data in the five SDSS optical bandpasses 
ugriz (jFukugita et al.l 11996). Only quasars that have 
nonzero flux in all five SDSS bandpasses after correcting 
for Galactic extinction have been kept in our catalog, 
bringing the sample size to 103,895. While the SDSS 
survey covers a large area of sky, it is limited to rela- 
tively bright quasars (i < 19.1 for z < 3 and i < 20.2 for 
z > 3). As such, in addition to these quasars, we included 
15,757 lower luminosity, optically-selected quasars taken 
from the Two Degre e Field QSO Reds hift survey (2QZ; 
bj < 20.85 for z < 3: lCroom et al|[200l . the Two Degree 
Field-SDS S LRG and Q SO survey (2SLAQ; g < 21.85 
for z < 3: IGroomet al 2009). and the AAT-UKIDSS- 
SDSS survey (AUS; i < 21.6 for < z < 5; Groom 
et al., in preparation). In the rest frame, we have over 
70,000 QSOs with log(i/[Hz]) ranging from -13.7-15.4 
(A -1200A-6/im) and at least 10,000 QSOs covering 
the range from ~13.3-15.6 (A ~750 A-15Aim). Ah SDSS 
magnitudes have been corrected for Galactic extinction 
according to lSchlegel et al.l (|1998[ ) with corrections to the 
extinction coefficients as given by iSchlaflv fc Finkbeineii 
(|2011f) . Table [1] presents the IR-UV photometry for the 
119,652 quasars in our sample. Signs of dust reddening 
are seen in 11,468 of these, specifically wi th a color ex- 
cess A{g -i) > 0.3 (|Richards et al.ll2003l ). These have 
been excluded from our analysis, bringing our final sam- 
ple size to 108,184. Figure [T] shows the sky coverage 
of our multi-wavelength data and Figure [2] shows the 
2500 A luminosity versus redshift distributions for each 
surve y. Note that , in th e so-called SDSS "Stripe 82" re- 
gion (jAnnis et al.ll201l[ ) along the Northern Equatorial 
Stripe, the SDSS imaging data reaches roughly 2 mag 
deeper than the main survey as a result of co-adding 
many epochs of data. 

2.1. Near- infrared 

Near-IR data in the JHK bandpasses are take n from 
the Two Micron AU Sky Survey f2MASS: ,Skrutskie et al.l 
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Fig. 1. — Footprints of the optically selected quasars in the matched surveys. Coverage of GALEX (blue), UKIDSS (teal), 2MASS 
(green), WISE (orange), and Spitzer (red) data are limited to the SDSS footprint (black) shown in the top-left panel. 



|2006[) . These include values that were matched to the 
2MASS All-Sky and "6x" point source catalogs using 
a matching radius of 2" with the 6x deeper catalog 
taking priority. In addition to these data, we include 
2cr extractions (i.e., "forced photometry") at the posi- 
tions of known SDSS quasars. These objects have suf- 
ficiently accurate optical positions that it is possible to 
perform aperture photometry at their expected locations 
in the near-IR imaging, despite being non-detections 
in 2MASS. Th e se 2a "detections" were cataloged by 
[Schneider et al.l (J201C0 : we include those objects with a 
signal-to-noise ratio (S/N) greater than 2. 

To supplement the 2MASS data, we have matched 
our catalog to the UKIRT (United Kingdom In- 
frared Telescope) Inf rared Deep Sky Survey (UKIDSS; 
[Lawrence et"aLl 12001 UK IDSS use s the YJHK fihe r 
svstem (jHewett et all [20060 : see also IPeth et~all (j2011[ ). 
In order to generate our combined optical and near-IR 
data sets, we begin by source matching samples of SDSS 
data with the UKIDSS LAS catalog using the Cross-ID 
form located on the Web site of the WFCAM Science 
Archive (WSAj]- In particular, we match against the 
UKIDSS DR5 LAS source table which contains the indi- 
vidual detections for a given object from each bandpass 

^ pittp : //surveys . roe ■ ac . uk : BOBO/wsa/crossID.f orm . j sp] 



merged into a single entry. We use a matching radius of 
0'.'7 and the nearest neighbor pairing option, accepting 
only the nearest object with a detection in at least one 
band as an acceptable match. 

Since UKIDSS is deeper than 2MASS, when a quasar 
has data in both surveys, we only use the UKIDSS data. 
Between UKIDSS and 2MASS, 77,864 of our sample 
quasars have data in the near-IR: 35,749 from UKIDSS 
data and 41,434 from 2MASS, including both 5a detec- 
tions and forced photometry. Figure |3] shows the dif- 
ference between the UKIDSS and 2MASS magnitudes 
versus the UKIDSS magnitudes for the 12,130 quasars 
the two surveys have in common. We have split each 
figure into three distinct groups: quasars taken from 
the 2MASS 5a catalog, quasars with 2a 2MASS extrac- 
tion values in more than one filter, and quasars with 2a 
2MASS extractions in only one filter. As the quasars in 
the third group have much larger scatter (up to 2.5 mag) 
that can result in a misestimation of the quasar's lumi- 
nosity by up to one order of magnitude, we have removed 
all quasars that fall into this group. This process leaves 
23,088 objects in our 2MASS sample and brings our total 
number of QSOs with near-IR data to 58,837. In Table 
[Hall 2MASS and UKIDSS fl uxes have be e n con verted to 
the AB magnitudes using Hew ett et al.l (|2006( ) conver- 
sions between AB and Vega magnitudes for the UKIDSS 
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Fig. 2. — Contour-scatter plot showing the 2500 A luminosity 
vs. redshift for each quasar color coded by survey. The contours 
indicate the linear density of scatter points on the plot with level 
drawn at 0.005, 0.5, and 0.95 of the normalized distributions. The 
histograms are normalized by the total number of quasars in each 
survey: 119,652 in SDSS, 42,043 in GALEX, 35,749 in UKIDSS, 
23,088 in 2MASS, 85,358 in WISE, and 1196 in Spitzer. The black 
dashed line represents the flux limit of the main SDSS survey before 
extending the sample to lower luminosities. 



system. 



2.2. Mid-Infrared 



To extend our SEDs into the mid-IR, we included 
matches to the WISE final data release (Wrigh t et alJ 
[2010). The W/5'i? bandpasses are generally referred to as 
Wl through 14^4 and have effective observed fi'ame wave- 
lengths of 3.36, 4.61, 11.82, and 22.13^m, respectively, 
for a typical quasar SED. The matching was performed 
by taking all non-contaminated WISE point sources 
within 2" of an SDSS quasar as a match. This matching 
radius maximizes the number of true objects matched 
while also minimizing the number of false matches. Fig- 
ure m shows the number of WISE matches found as a 
function of separation distance from SDSS quasars. To 
ensure only the best matches were included, we also re- 
quired that all matched quasars have S/N > 10 in both 
Wl and W2. In total, we have 85,358 matches to WISE; 
we estimate the false matching rate to be 1% based on 
the number of matches obtained after shifting all WISE 
data by 200" in declination (red line in Figure S]) . Since 
the WISE photometry was calibrated against blue stars 
it tends to overestim ate the flux of red s ources by about 
10% in the WA band (| Wright et al.ll2010f ). this correction 
has been applied to all quasars in our sample. 

When available, Spitzer IRAC data were also in- 
cluded. We specifically included dat a from t he Extra- 
galactic First Look Survey (XFLS; iLacv et al. 2005}; 
Spitzer Deep-Wide Field Survey (SDWFS; Aslib v et atl 
mm : SWIRE DR2 (Lonsdale ct al. 2003 ), including the 
ELAIS-Nl, -N2, and Lockman Hole fields; S-COSMOS 
(jSanders et al.l [20071 ). and our own extraction of high 
redshift QSOs in stripe 82 (S82IIIZ, program number 
60139). The fluxes for known SDSS quasars from the 
Stripe 82 program are reported here for the first time in 



Tahle\B [Spitzer source labeled as S82HIZ). The IRAC 
bandpasses (chl-ch4) have effective wavelengths of 3.52, 
4.46, 5.67, and 7.7 ^m for a typical quasar SED. There 
are a total of 1196 matches to Spitzer. For a small num- 
ber of quasars ('^ 200) both Spitzer and WISE data are 
available; in these cases, we use both data sets. In Ta- 
ble [1] all Spitzer and WISE fluxes are reported in AB 
magnitudes. 

2.3. Ultraviolet 

To extend the SEDs into the UV, we use G alaxy 
Evolution Explorer (GALEX; iMartin et all 120051 ) data 
when available. The effective wavelengths for the near- 
UV (NUV) and the FUV bandpasses are 2267A and 
1516A, respecti v ely. M ost of our matches are taken from 
iBudavari et"al[ (J2009I ) who matched GALEX-GR6 to 
SDSS-DR7. For this study, we have taken only the most 
secure matches; that is, only one SDSS quasar matched 
to only one GALEX point source. This algorithm gives 
14,302 matches and leaves out 53,782 quasars that have 
multiple SDSS sources matching to the same GALEX 
source. In addition to these data, we inc lude forced 
point - sprea d function (PSF) photometry (jBovv et al.l 
120111 I2012f) on the GALEX images in the positions of 
the SDSS point sources, which adds 61,490 detections. 
As with 2MASS, we onl y kept the extracted data with 
S/N > 2, and as with the lBudavari et all ()2009n matches, 
we limit the sample to those with no source confusion; 
this reduces the number of extracted quasars to 27,744. 
Combining both sets, the total number of matches is 
42,046. All GALEX photometry has been corrected for 
Galactic extinction, assuming v4nuv — 8.741 x E{B — V) 
and yJFUv = 8.24 x E{B - V) ~ 0.67 x E{B - Vf 
(iWvder et al.ll2007l) The^ (B-y) values are taken from 
the lSchneider et ahl (J2010I ) catalog. All the matched data 
are given in Table [U 

2.4. X-Ray 

To construct full SEDs, we must extend our data 
into the X-ray regime. Due to comparatively limited 
sky coverage of sensitive X-ray observations, the num- 
ber of our quasars with X-ray data is quite small com- 
pared to the size of the samples in the optical and IR. 
For exam ple, there are only 277 matches in the ChaMP 
data set (jGreen et al.ll2009h . So instead we take advan- 
tage of the careful work on mean X-ray p roperties as a 
function of UV luminosity as compiled by iSteffen et alJ 
(2006) ■ Speciflcally, we determine the X-ray flux of each 
quasar using the Luy-Lx relation that parameterizes 
the correlation between the 2500A and 2keV luminos- 
ity of quasars. We flnd the 2500 A luminosity by ex- 
trapolatio n from the close st fllter using an ctopt = —0.44 
(Vandcn B erk et al.l[2001h . The 2keV luminosity is es- 
timated using Equation ^ and our 2500 A luminosities. 
An X-ray energy spectral index Ofa, = — 1 (photon in- 
dex r = — 2) was assumed between 0.2 keV and lOkcV 
(e.g. iGeorge et al.l [2000l ) . In this way, we can esti- 
mate the X-ray part of the SED for all of our sources 
rather than using just the small fraction of objects with 
robust X-ray detections. This process ignores the cor- 
relations betw e en Qq x (and a^v) and a^ as discussed in 
iKruczek et al.l (|2011f ); however, these trends are small as 
compared to the overall trend with luminosity. 
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Fig. 3. — Contour-scatter plot showing the difference between the 2MASS and UKIDSS magnitudes vs. the UKIDSS magnitudes for the 
J and K filters. The blue contours/dots show the catalog 2MASS values, green shows the 2MASS 2cr extraction values with a detection 
in only one band, and red shows the rest of the extracted values. The blue dashed line shows the 2MASS point source catalog's limiting 
magnitude. The two surveys contain 12,130 quasars in common. Green points were excluded as the errors are quite large (see text). Similar 
results were found for the H filter. 

through the fiher curves as discussed below. For the host 
galaxy correction, we subtract a model template. We do 
no t consider beaming specifically, but refer the reader 
to iRunnoe et al.l (|2012[ ) for a discussion of how our re- 
sults would chang e under the assumption of non- isotropic 
emission; see also lNemmen fc BrothertonI ()201ClD . 

3.1. Lyman Forest and Limit 

A particular challenge is to determine the SED in the 
extreme-UV (EUV) (Arost < 1216A) part of the spectrum 
where intergalactic hydrogen causes significant attenua- 
tion of the quasar signal (Lya forest; Lvnds'1971f). To ac- 
count for this attenuation, the redshift-depcndent effec- 
tive optical depth. Teg (z, A), of Meiksin (2006) was used. 
This optical depth models the average attenuation of a 
source assuming Poisson-distributed intergalactic hydro- 
gen along the line of sight, out to a given redshift. This 
optical depth is split into three parts: the contribution 
due to resonant scattering by Lyman transitions, systems 
with optically thin Lyman edges, and Lyman Limit Sys- 
tems (LLS). Figure [5] shows this model for redshifts of 1, 
2, and 3 for both LLS and non-LLS. 

Because our investigation uses photometry and not 
spectra, we do not know whether a LLS is present. Fur- 
thermore, we only have spectral coverage of the LLS re- 
gion for a fraction of our sources. Thus, to be as con- 
servative as possible, we assume that there is no LLS 
present; that way we only make the minimum correc- 
tion needed for each SED. However, it is important to 
note that the SDSS quasar selection algorithm has been 
shown to preferentially sele ct quasars possessing a LLS 
in the range 3 < 2 < 3.5 (Worseck fc Prochaska 2011) . 

For the continuum {Fy) we have assumed a power- 
law of the form i'~°'^'^ (iVanden Be rk et al.l l200lD . consis- 
tent with the results of iScott ~ al. (2004). Continuum 
weighting is necessary since we are dealing with broad- 
band photometry and not spectroscopy. With high res- 
olution spectroscopic data the correction is exact as the 
iMeiksinI (|2006f) corrections themselves are independent of 
the SED. However, for broadband photometry, where the 
bandpass can overlap features in the r distribution, the 
convolution of the SED with the filter response changes 
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Fig. 4. — Angular separation between known quasars and nearby 
WISE sources (blue) and the known quasars with declinations 
shifted by 200" (red). The green curve shows the best fit expected 
distribution (Rayleigh distribution for small separations and linear 
growth for larger separations). The vertical magenta dot-dash line 
indicates the matching radius used. As seen in the inset, the blue 
histogram does not follow the expected distribution; this is due, in 
part, to the large beam size of WISE (FWHM=6"). 

As a check we have compared ou r X-ray extrapolations 
to ChaMP data from lGreenet al.l (|2009l see Section g]). 
In all cases we found the X-ray data to fall within the 
1.5a limits of the ttox extrapolation. 

3. CORRECTIONS 

When studying quasar SEDs, we are interested in the 
true continuum level of the radiated light. The con- 
tinuum may be contaminated by, for example, spec- 
tral emission lines, absorption by intergalactic hydro- 
gen clouds, host galaxy contamination, and beaming 
effects. Here we address each of these features. For 
the hydrogen absorption and broad emission lines, we 
determine a magnitude correction by folding a model 
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Fig. 5. — IMeiksiiJ l|200g ) model for Lyman series extinction as 
a function of observed wavelength for z = 1 (blue), z = 2 (green), 
and 2 = 3 (red). The solid lines show the extinction assuming a 
LLS (at the corresponding redshifts) along the line of sight, and the 
dashed lines show the extinction without a LLS. The jagged edges 
are a result of absorption from the Lyman series and the drop-off is 
due to absorption from the IGM and any hydrogen clouds along the 
line of sight. The transmission rises with decreasing wavelength for 
each feature because the shorter wavelengths sample lower redshifts 
and the universe is more ionized at lower redshifts. 



the effective wavelength of the bandpass. 

To apply the Lyman forest correction, we convolve a 
continuum with the IMeiksinI ([2006) optical depth and 
each filter over < z < 6 in Az steps of 0.01 and use lin- 
ear interpolati on to precisely match redshifts. Following 
IMeiksinI ()2006f ) , the correction is given in magnitudes as 
a transmission weighted average of r: 



K 



IGM 



= -2.5 X log 



/ \SxFxe-~'-''''d\ 
/ XSxFxdX 



(3) 



where Sx is a transmission filter curve, -Fa is the contin- 
uum, and Teff is the redshift-dependent effective optical 
depth. This value is calculated independently for each 
filter and makes m{T = 0) = mir) — iiTiGM such that 
all quasars are made brighter in the Lyman absorption 
regions as a result of the correction. 

3.2. Emission Lines 

The presence of emission lines in a photometric 
bandpass affec t s quas ar magnitudes as illustrated by 
IRichards et al.l ()2006bl e.g., see their Figures 8 and 17). 
Thus, to make a fair comparison between sources at all 
redshifts, we endeavor to remove first-order effects of the 
emission line contributions to the measured magnitudes. 
Our emission-line template takes the position, width, and 
equivalent width (EW) of th e 13 strongest spectral line s 
(as labeled in Figure [S]) from lVanden Berk et al.l (|2001h . 
We make no attempt to correct for the small blue bump 
(i.e. Balmer continuum or Fe II emission), but recog- 
nize that those features can have a significant impact; 
indeed their residuals can be seen in our mean SEDs. 
Figure |6] shows the mock spectrum used in our emission- 
line corrections; it includes a double power-law contin- 
uum {uy = —0.46 for Arcst < 4600 A and a^, — —1.58 for 
Arest > 4600 A) and 13 emission lines. 
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Fig. 6. — IVanden Berk et al.l 11200 II ) composite spectra (black) 
with our model spectrum with a LLS (red) at a redshift of 3.75. 
Since the SDSS spectra cover 3800—9200 A, the contributors to the 
points at Ajost < 1216 A are high-^ quasars and are consistent with 
the model used. Since we model the main emission line features, on 
average, they should not appear in the SEDs. However, unmodeled 
residuals in Fe II and the Balmer continuum (small blue bump) can 
clearly be seen in the SEDs on either side of the Mg II emission line 
at 2800 A. We use a double power law continuum in this figure since 
the composite spectra do not correct for host galaxy contributions 
on the red end. Model spectra used to correct our SEDs are placed 
on a single power law as shown in Figure [Tj 



The correction is computed as 



Xcm = -2.5 X log 



J XSxFxicM)dX \ 
JXSxFx{c)dX )' 



(4) 



where Sx is a transmission filter curve, Fx (c) is the con- 
tinuum, and Fx(ck.l) is the continuum with spectral 
lines. This value is calculated independently for each 
filter and makes m{c) = m{cSzl) — i^cm, such that all 
quasars are made fainter by the removal of the emission 
line contribution. 

We further use this template to illustrate the correc- 
tion for hydrogen absorption (assuming the quasar is at 
z = 3.75 and there is a LLS). Overall, we find reason- 
able agree ment (by eye) with the rn ean composite spec- 
trum from lVanden Berk et al.l (|2001l ) in terms of the con- 
tinuum, broad emission lines, and intergalactic medium 
(IGM) attenuation. 

In Figure [71 we illustrate the effect of these emission- 
line and IGM corrections using a mock spectrum of a 
z = 3 quasar (without a LLS). Colored points indicate 
the observed instrumental magnitudes and the intrinsic 
continuum magnitudes after applying the emission line 
and IGM ivT-corrections. 

3.3. Host Galaxy 

It is also necessary to correct the data for host galaxy 
contamination. Since we lack the data to measure this 
directly, we must estimate the host galaxy contribution. 
To do this, we combine two differ ent models: the rela- 
tion outlined in iShen et all ([20111) for t he higher lumi 



nosity quasars, and the relationship from Richards et al. 



(,2006a,) for the lower luminosity quasars. The .Shcn et al. 
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Fig. 7. — Observed frame mock spectrum (solid line) and Qj, = 
—0.44 continuum (dotted line) for a 2 = 3 quasar. The closed dia- 
monds bifurcated by solid lines indicate the "measured" instrumen- 
tal magnitudes, while open diamonds bifurcated by dashed lines 
indicated the corrected values (assuming an underlying power-law 
SED). The width of the lines indicate the width of t he b andpasses. 
Here we have correcte d for IGM attenuation (SectionlsTTJ and emis- 
sion line flux (Section [32J. 



()2011| ) relationship is 



5iooA,host ^ 0.8052 - 1.5502a; + 0.91212;^ 

5100A,QSO 

where 



0.1577a;^ (5) 



log 
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■44 



and 

< a; < 1.053, 

which sets the relative scaling at 5100A. To extend this 
expression t o all wavelengths, we use th e ellipt ical galaxy 
template of iFioc fc Rocca-Volmeraiigi (|1997D scaled to 
5100A. When x > 1.053, the quasar completely outshines 
the host galaxy and no correction needs to be applied. 
We note that this relationship was found using quasars 
that have log (A °L » ) > 44, and cannot, in 

° ^ 5100A SlOOA.totaF — ' ' 

general, be extrapolated to lower luminosities. 

W e found t ha t, fo r lower luminosity quasars, 
the iShen et all ()201lD equation overestimates the 
host galaxy; we ins tead use the relationship from 
[Richards et al.l ()200 6a') which, was adapted from the re- 
lationship used byLVandcn Berk ct al. (2006): 



-H2.887-log(LBoi/iEdd) 



(6) 



where iog{L^^^^^^^^^) = \og 



G156A,total 



erg s ^Hz' 



■ 1 



6156A, host 
erg s~lHz~^ 



This system of equations is then solv ed numerically 
for the host galaxy luminosity. As in [Richards et alJ 
(|2006af) we take LBoi/^Edd to be unity, since this pro- 
vides the minimum correction needed. This sets the rel- 
ative scaling at 6156A and, a s before, we use the ellipti- 
cal galaxy template of lFioc fc Rocca-Volmcrangc, (, 199 7,1 



to extend this to all wavelengths. To ensure a smooth 
transition between these two methods, we have chosen 
the crossover luminosity to be the point at which both 
methods agree, log (A °L ° ) — 44.75. As host 

^ ' o \ 5100A 5100A, total' 

galaxy subtraction can have a large impact on the SED 
near 1 /xm and it is not a well-esta blished procedure, w e 
not e that other methods inc lude iCroom et al.l ()2002D . 
and lMaddox fc Hewettl (|2006f l. 

It was brou ght to our attention th at the host galaxy 
model used in iRichards et al] (|2006a| ) was incorrect (K. 
Leighly, private communication 2012). The model was 
converted from F\ to F^, but was not converted to vFi, 
before subtracting it from the SEDs. This does not effect 
the SED at 5100A, where the SEDs are normalized, but it 
is systematically wrong on either side of this wavelength. 
Herein we correct that error. 

3.4. Gap Repair 

Since all the quasars in our sample have been selected 
from SDSS photometric data, they all have ugriz mea- 
surements, but are not guaranteed to have measurements 
in the other bandpasses that we utilize. To address this 
issue we will "gap repa i r" all m issing data in a way sim- 
ilar to IRichards et al.l (|2006aD . Specifically, we replace 
missing values with those determined by normalizing an 
interpolation/extrapolation of the continuum in the next 
nearest bandpass for which we have data. Sometimes this 
will be a previously constructed mean SED; other times 
it will be a functional form. To estimate errors for the 
gap filled points, we have fit up to third-order polyno- 
mials to plots of the magnitude errors, log(crm), versus 
magnitudes, m, for each of the filters. These functions 
were then used to estimate Um for the gap filled value for 
m. 

This procedure works well where we can interpolate be- 
tween filters (e.g., at wavelengths longward of the Lyman 
limit). Beyond the Lyman limit, we are no longer inter- 
polating, but extrapolating, which is a trickier process. 
Using limiting magnitudes does not solve the problem as 
redshift effects mean that even if all quasars had coverage 
from all of the filters, the same rest-frame wavelengths 
are not observed for both high- and low-redshift quasars. 

On the long wavelength end, we always extrapo- 
late to longer wavelen gths using the mean SED from 
IRichards et al.l ()2006aD . normalized to the nearest mea- 
sured bandpass. On the short wavelength end, the gap 
repair process depends on whether we are making a mean 
SED or attempting to reconstruct the SED for an individ- 
ual quasar (e.g., to determine its bolometric luminosity). 
See below for details for specific cases. 

Unless otherwise stated, the X-ray part of the SED 
is de termined usin g Equation ([2]), the iuv~aox relation 
from lSteffen et ahl 12006). and linearly interpolating be- 
tween the X-ray and the bluest data point after gap re- 
pair. The errors in this region are dete rmined from the 
uncertainties of the iSteffen et al.l ()2006[ ) relation. 

4. MEAN SEDS 

4.1. Overall Mean SED 

To determine the mean quasar SED, we converted all 
the flux densities for each quasar to luminosity densi- 
ties and shifted each broadband observation to the rest 
frame. The data were then placed onto a grid with 



points separated by 0.02 in log (ly); kriging was then used 
to align the broadband luminosities to our grid points. 
Kriging is a nonparametric interpolation method that 
predicts value s and errors for r e gions between observed 
data points. IRvbicki fc Press! ( 19921 ) were among the 
first to present this technique in the astronomical lit- 
erature (under the name Wiener filtering). Since then 
it has been mainly u sed to estimate light curves (e.g., 
iKozlowski et al.ll2010[ ). 

Kriging estimates the correlation between data points 
as a function of separation, then uses this correlation 
to interpolate between the data points. We used the 
R package gstat to perform variance- weighted kriging 
using an exponential variogram model. Once the data 
are rebinned in the rest-frame, the arithmetic mean and 
standard deviation are taken at each grid point. Because 
kriging estimates variance values at each of the new grid 
points, we are able to estimate errors on our bolometric 
luminosities and the derived bolometric corrections. 

Figure |8] shows the "raw" data points (without the 
corrections described in Section [3]) and mean for all the 
quasars in our sample. On the red end, the SED drops 
due to a lack of data and on the blue end, shortward of 
the Lya line, there is a clear drop due to intergalactic ex- 
tinction. The flattening in the mean shortward of the Ly- 
man limit is due to redshift effects; low-redshift sources 
have no rest-frame measurements at this frequency, bi- 
asing the data toward higher luminosity. 

To find the corrected mean SED, we first construct a 
"gap fiUing" SED by looking at the mean SED of - 2100 
quasars with full wavelength coverage; i.e., they have 
at least four data points in mid-IR, at least three data 
points in the near-IR, and full coverage with GALEX 
in the UV. Redshift effects cause the SED coverage in 
the rest frame to drop off sharply around the Lyman se- 
ries. To avoid the mean being biased towards the higher 
luminosity SEDs in this region we have truncated this 
gap filling mean SED where the total number of quasars 
drops (around 1216 A) and use the mean UV luminosity, 
{L °), to find the X-ray luminosity using Equation 

([2]). Figure m shows the resulting mean SED that we will 
use for gap filling individual quasar SEDs that lack full 
coverage in all of the bandpasses considered herein. 

To find the mean SED for the entire sample, we apply 
the corrections described in Section [S] but only gap re- 
pairing filters with wavelengths A > 912A, so that way 
we avoid gap filling in a region that is not well sampled 
with our gap-filling SED. We then truncate the mean at 
912 A and use the mean UV luminosity, {L »), to find 

the X-ray luminosity using Equation ([2]). We do this in- 
stead of using actual X-ray data since the re are too few 
X-ray detections and Steff en et al.l (p006) have already 
done the careful work of extracting the relationship be- 
tween the UV and the X-ray parts of the SED. 

Both radio-loud and radio-quiet quasars are included. 
The mean SED of 108,184 quasars is given in tabular 
form in Table [2] Figure [TU] shows the resulting mea n with 
t he Ri chards et al. (2006a) mean, the Vanden Ber k et al.l 
(|2001t ) power law, a typical host galaxy, and the ChaMP 
X-ray data included for compari son. The difference be- 
tween our mean SED and that of lRichards et al.l (|2006a[ l 
in the 9 12-12 16 A region is due to our more careful cor- 
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Fig. 8. — Top: mean uncorrected SED (solid line), data points for 
all quasars, and la level (dash-dotted line) for the raw sample (no 
corrections applied, see Section [Sj. The steep drop-off on the red 
end is caused by the lack of high luminosity, low redshift quasars 
in our sample; this is caused by both quasar evolution and the fact 
that at lower redshifts we sample a smaller volume of the universe. 
The white wedges between the filters are caused by the limiting 
magnitudes for each filter. The drop-off just shortward the Lya 
line is due to IGM attenuation, but the flattening shortward of 
the Lyman limit is due to the lack of low-luminosity, high-redshift 
quasars, caused by the detection limits of the surveys. Bottom: 
the number of SEDs averaged at each frequency. The two vertical 
magenta lines show the positions of the Lyo line at 1216 A and the 
Lyman limit at 912 A. 

rections of systematic effects in this region as described 
in Section [31 The abrupt change in the slope of the SED 
at '^ 1100 A is not expected to be real, but rather repre- 
sents our lack of knowledge of this region of the SED. We 
are largely limited to simply connecting the two better 
known regions of the NUV and soft X-ray with a power 
law. We will further discuss the range of possible FUV 
continua in Section [6l 

4.2. Sub-sampled Mean SEDs 

While the overall mean quasar SED is a useful tool, 
examining how the SED changes as a function of various 
quasar parameters may shed light on th e physical pro- 
cesses of the central engine. For example. [Richards et al.l 
(2006a) found that bolometric corrections (see Section[5]) 
differed by as much as a factor of two in the extremes of 
quasar types, but with only 259 objects [Richards et al.l 
(|2006al) did not have the data necessary to comment 
on what ph ysics is behin d the r ange of bolometric cor- 
rections. M arconi et al.l (|2004[) provide a luminosity- 
dependent bolometric correction, but this is largely de- 
pendent on the J^uv~Q^o x relationship and is built into the 
[Richards et al.l ([2006aD bolometric corrections. What we 
seek (with the aid of a substantially larger sample) is a 
more physical understanding of the SED differences and 
the resulting changes in bolometric correction. 

In an attempt to better understand the physics that 
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TABLE 2 
Mean Quasar SEDs 



All 



CAll 



Low Lum. 



Mid Lum. 



High Lum. 



C IV Line 
Zone 1 



C IV Line 
Zone 2 



UV Bump 



Scott UV 



Casebeer UV 



13.0 


45.201 


0.52 


44.844 


0.24 


45.272 


0.2 


45.632 


0.25 


45.434 


0.31 


45.831 


0.33 


45.213 


0.55 


45.248 


0.55 


45.248 


0.55 


13.02 


45.221 


0.52 


44.864 


0.24 


45.292 


0.2 


45.652 


0.25 


45.454 


0.31 


45.851 


0.33 


45.233 


0.55 


45.268 


0.55 


45.268 


0.55 


13.04 


45.231 


0.52 


44.874 


0.24 


45.302 


0.2 


45.662 


0.25 


45.464 


0.31 


45.861 


0.33 


45.243 


0.55 


45.278 


0.55 


45.278 


0.55 


13.06 


45.241 


0.52 


44.884 


0.24 


45.312 


0.2 


45.672 


0.25 


45.474 


0.31 


45.871 


0.33 


45.253 


0.55 


45.288 


0.55 


45.288 


0.55 


13.08 


45.251 


0.52 


44.894 


0.24 


45.322 


0.2 


45.682 


0.25 


45.484 


0.31 


45.881 


0.33 


45.263 


0.55 


45.298 


0.55 


45.298 


0.55 


13.1 


45.271 


0.52 


44.914 


0.24 


45.342 


0.2 


45.702 


0.25 


45.504 


0.31 


45.901 


0.33 


45.283 


0.55 


45.318 


0.55 


45.318 


0.55 


13.12 


45.281 


0.52 


44.924 


0.24 


45.352 


0.2 


45.712 


0.25 


45.514 


0.31 


45.911 


0.33 


45.293 


0.55 


45.328 


0.55 


45.328 


0.55 


13.14 


45.291 


0.52 


44.934 


0.24 


45.362 


0.2 


45.722 


0.25 


45.524 


0.31 


45.921 


0.33 


45.303 


0.55 


45.338 


0.55 


45.338 


0.55 


13.16 


45.301 


0.52 


44.944 


0.24 


45.372 


0.2 


45.732 


0.25 


45.534 


0.31 


45.931 


0.33 


45.313 


0.55 


45.348 


0.55 


45.348 


0.55 


13.18 


45.301 


0.52 


44.944 


0.24 


45.372 


0.2 


45.732 


0.25 


45.534 


0.31 


45.931 


0.33 


45.313 


0.55 


45.348 


0.55 


45.348 


0.55 



Note. — All of the SEDs are taken to have a^ — —1 above 0.2 keV. Units are log(erg s ). This table is published in its entirety in the electronic edition of 
the online journal. A portion is shown here for guidance regarding its form and content. 
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Fig. 9. — Top: mean "gap filling" SED (solid line), based on 
~2100 quasars with full broadband wavelength coverage. To avoid 
a bias towards high luminosity at the blue end where the number of 
quasars drops off, the mean is truncated at 1216 A, and connected 
with the X-ray (see text). Bottom: the number of SEDs averaged 
at each frequency. The dotted line indicates the data that is not 
used due to truncation (see text). The two vertical magenta lines 
show the positions of the Lya line at 1216 A and the Lyman limit 
at 912 A. 

lead to differences in SEDs (and bolometric corrections), 
we consider two parameters herein, specifically the UV 
luminosity and C IV emission line properties. While it 
is also important to consider SEDs as a function of mass 
and accretion rate, those quantities are not directly mea- 
surable; we will leave that analysis to future work (C. 
M. Krawczyk et al. 2013 in preparation). Luminosity- 
dependent SEDs (Section 14.2. ip are of interest because 
of the known Luv~Q;ox relationship and strong depen- 
dence of accretion disk wind physics on the UV to X- 
ray flux ratio (e.g. iProga et al.ll2000l ). Examining the 
mean SED as a function of C IV emission line properties 
fSection l4.2.2l) is interesting because that line may be an 
indicator of the true SED as it serves as a diagnostic of 
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Fig. 10. — Top: mean corrected SED (black) and data points 
for 108,184 quasars with gap filling, edge filling, host galaxy 
removal, and spectral corrections (see text). The mean SED 
from Richards et al. (2006a, cyan), an a^ = —0.44 power law 
(Vanden Berk et al. 2001, dashed grey), a typical elliptical host 
galaxy (green), and ChaMP X-ray data (Green et al. 2009, black 
crosses) arc included for comparison. The mean has been truncated 
at the Lyman limit to avoid being biased to high-luminosities when 
the number of quasars drops off (see text). See legend of Figure [51 
for color coding of scatter points. Bottom: the number of SEDs 
averaged at each frequency. The dotted line indicates where the 
mean SED is truncated and connected with the X-ray. The two 
vertical magenta lines show the positions of the Lyo line at 1216 A 
and the Lyman limit at 912 A. 

which two com ponents dominate the broad emission lin e 
region (BELR; iRichards et all [20TTI : iWang et al.llMll ). 
In fact, UV emission lines like C IV, with ionization po- 
tentials in the EUV part of the spectrum, may even be 
an indicator of the unseen EUV part of the SED (see 
Section [6]). 

4.2.1. Luminosity-dependent Mean 

The well-known nonlinear relationship between the UV 
and X-ray luminosities means that the SED (and thus 
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Fig. 11. — Cumulative histogram of loe(i^Ly) I • for all 

non-reddened quasars. Vertical lines show luminosity cuts that 
result in three equally populated bins. The inset shows the redshift 
distribution for these three bins, with the overlap shaded in. 



the bolometric correction) must be a function of lumi- 
nosity. As a first test, we have spht the sample into 
three equally populated luminosity bins, each contain- 
ing 36,061 quasars. Figure [TT] shows the cumulative his- 
togram of log {vL^) I ° for all quasars, with the 

vertical lines showing the luminosity cuts. In order to 
separate the changes that potentially arise from evolution 
from those that are luminosity-dependent, we marginal- 
ize over redshift by selecting sub-samples from each bin 
that each have the same redshift distribution. The inset 
of FigurefTTJshows the redshift distribution for each of the 
bins; the shaded region shows the redshift distribution of 
the sub-samples. This process leaves ~10,000 quasars in 
ea ch luminosity b i n as co mpared to the 259 total objects 
in [Richards et all (|2006a| ) across all luminosity bins. 

We calculated the luminosity-dependent mean SEDs in 
an iterative way. As a first step we follow the same steps 
as the overall mean SED: we used our gap-filling SED to 
interpolate/extrapolate over gaps in the photometry with 
A > 912 A and truncated the mean SEDs at 912A. These 
mean SEDs were then connected to the X-ray points de- 
termined using the mean UV luminosity for each sample. 
The resulting mean SEDs were then used as the new gap- 
filling SEDs for each luminosity bin. For all remaining 
iterations all missing photometry was gap filled, and the 
resulting mean SEDs were truncated at 800 A (where the 
sampled number starts to fall) before connection to the 
X-ray. Figure [12] shows the mean SEDs and data for each 
of the luminosity bins after 10 iterations. The iteration 
process described above ensures that our mean SED fol- 
lows the quasars with data in each luminosity bin instead 
of the initial gap filling model. The redshift marginaliza- 
tion process means that the mean redshift of all three 
luminosity samples is 2 ~ 1.5. The mean SEDs are given 
in tabular form in Table [5] 

Figure [T^ a) shows the mid-IR region of the three 
SEDs when they are normalized at 1.3/im and Fig- 
ure [TSjb) shows the optical (A < 3000A) and UV region 



(A > 3OOOA) when normalized at 1450A. The thick hues 
indicate where a Welch's t test shows the mean SEDs 
have less than a 1% chance of being the same {p <0.01). 
The high-luminosity SED has more hot dust emission 
(at 2-4/im), a small Balmer continuum (A < 3000A), a 
harder (bluer) optical spectrum, and a softer (redder) 
UV spectrum. The low-luminosity quasars would have 
to have host galaxies ~8 times more luminous than as- 
sumed in order for the optical slopes at A < 3000A to 
agree between the low-luminosity and high-luminosity 
SEDs (thus the redder continuum is likely to be intrin- 
sic and not due to host-galaxy contamination). The UV 
region is sensitive to the Lyman series corrections (Sec- 
tion 13. 2p . but these corrections are not luminosity de- 
pendent. The Welch's t test shows the prominent 10 /im 
silicate bumps in the high-luminosity objects are statis- 
tically si gnificant. The IR bum ps are discussed in more 
detail in lGallagher et all (|2007f ). 

One of our goals in considering the luminosity- 
dependent mean is to investigate the shape of the SEDs 
over 500 < Ajost (A) < 1200. Given the corrections nec- 
essary in this region (Section [3^ . the shape here is nec- 
essarily uncertain; however, the data supports a harder 
slope with lower luminosity (see Figure [TSTblV This is 
consistent with the results of lScott et al.l (|2004[ ) and sug- 
gests that the X-ray fiux may not be the only part of the 
SED that has a nonlinear relationship to L 



2500 A 



This 

behavior is also se en in the spectral Principal Component 
Analysis (PGA) of lYip et all (|2004D who found that their 
third eigenvector shows an anti-correlation of the con- 
tinua on either side of Lya. In short, quasars that are 
bluer longward of Lya have softer (redder) continua in 
the Lya forest region. This can have important conse- 
quences for bolometric corrections; see Section [5l 

4.2.2. C IV-dependent Mean 

While the Baldwin Effect (|Baldwinlll977l ) reveals that 
there is a relationship between the stren gth of UV emis- 
sion li nes and the continuum luminosity, iRichards et al.l 
()2011|) have argued that the G IV (and other UV emis- 
sion lines) properties are better diagnostics of the shape 
of the SED than its absolute scaling. If that is the 
case, SEDs made as a function of emis sion line prop- 
erties such as C IV blueshift and E^ y (IRichards et ah 



2011 1 ) or PGA (e.g., "Eigenvecto r 1" : iBoroson fc Green 



1992t iBrotherton fc Francis"1999^ may reveal interesting 



differences. Any such differences would have important 
implications for bolometric corrections of individual ob- 
jects. 

Thus, in addition to luminosity sub-samples, we have 
also divided the data as a function o f G IV blueshift and 
rest-frame EW ([Richards et al.ll201lD . To see ff the SED 
shape depends on these properties, we have taken two 
zones that are represen tative of extrema in disk-wind 
structures according to IRichards et al.l ([20111 see also 
[Wang et al.l ([20TT[ )) : Zone 1, blueshifts < 600 km s'^ 
and EW > 50 A (i.e., with "disk-dominated" BELRs) 
and Zone 2, blueshifts > 1200 km s'^ and EW < 32 A 
(i.e., with "wind-dominated" BELRs). These cuts are 
chosen to have roughly the same number of objects in 
each sample. In Zone 1 there are 5736 quasars with 
a mean losivLj^) I °= 45.9, and Zone 2 contains 

°^ '^' 'A=2500A ' 
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Fig. 12.— Mean SEDs for low {\og{uL„) | 
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(d) Luminosity dependent SEDs 



A— 45.41) luminosity quasars (a), mid (45.41 < log(t'Li/) | ,„;< 45.85) (b), and 

high (45.85 < log (i/Li/) | „„r) (c). Panel (d) shows the three curves overplotted. See Figure [8] for an explanation of the scatter points in 

A — ^DUUA. 

panels (a)-(c) and the bottom histogram of panel (d). The high luminosity SED has more hot dust emission (at 2-4/im), a smaller Balmer 
continuum (A < 3000A), a harder (bluer) optical spectrum (A < 3000A), and a softer (redder) UV spectrum (A > 3000A); see Figure [131' a) 
and (b). 



5713 quasars with a mean log [vL^) \ ° = 46.2. Fig- 

ure [14] shows the data and mean SEDs for the quasars 
in each zone. For these mean SEDs, we used the mean 
mid-luminosity SED for the initial gap repair to all the 
filters and truncated the mean at 8OOA before connect- 
ing it with the X-ray. We then took the resulting mean 
SED and used that to gap repair in the same way; this 
was repeated 10 times. We note that there is a difference 
in mean luminosity between Zones 1 and 2; however, we 
have not gap-filled with the luminosity-dependent SEDs 
in order to isolate any residual differences. The mean 
SEDs are given in tabular form in Table [2] 

Our primary interest in these SEDs lies in the UV part 
of the spectrum; however, it seems that there is rela- 
tively little that can be learned about that region from 



this sample. Seeing C IV requires that the sample is re- 
stricted to 2; > 1.6, which has a number of consequences. 
In particular, it limits the sample to the highest luminosi- 
ties and the reds hift limit introduces a bias toward LLS 
(jWorseck fc Prochaskall2011[ ). which may be the cause 
of the steep fall-off beyond Lya. However we are able 
to see differences in the optical/ UV continuum consis- 
tent with eigenvector 1 analyses (jBrotherton fc Francis! 
Il999f l. In particular. Figure [tsT c) shows that the Zone 1 
mean SED has a larger Balmer continuum (A < 3000A) 
and stronger Lya emission (the difference in slopes is 
consistent with being due to stronger C IV and Lya in 
Zone 1 mean SED relative to Zone 2 since we have only 
corrected for the mean emission line contribution). We 
will further consider differences in the bolometric cor- 
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Fig. 13. — Normalized SEDs with the y-axis showing relative luminosity. The thick lines indicate where the probability of the mean 
SEDs being the same is less than 1% and the vertical dashed lines indicate the normalization frequency, (a) Luminosity-dependent SEDs 
normalized at 1.3 /xm; high-luminosity quasars appear to have more hot dust emission, (b) Luminosity-dependent SEDs normalized at 
1450 A; high-luminosity quasars have softer EUV continua. (c) C IV-dependent SEDs normalized at 2027 A. Zone 1 ( disk-dominated) 
quasars have more Balmer continuum and stronger Lya emission then Zone 2 (wind-dominated) quasars. See Figure [141 (d) Bolometric 
correction-dependent SEDs normalized at 2500 A, showing the effect that most of the range in BC comes from the range of UV and EUV 
continua (and dust). 



rection and wliat can be learned about the sliape of the 
SED as a function of the UV emission hne properties in 
Section [5l further investigations of the UV continuum of 
quasars as a function of C IV emission hne properties are 
certainly warranted. 

5. BOLOMETRIC CORRECTIONS 

One of the main goals of this study is to characterize 
the bolometric luminosities of quasars. The bolometric 
luminosity is the integrated area under the SED curve 
and can be calculated as: 



L 



bol 



L,,dv 



hi{\{i)vL^d\og{v) (7) 



where Figures [5] and [TO] plot the quantity vL^. Nor- 
mally, we do not have an accurate measurement of the 
full SED for a quasar, but we can estimate the bolometric 
luminosity by assuming some SED template and using a 
known monochromatic luminosity, Ly. This "bolometric 
correction" (BC) is given as: 



BC, 



L 



bol 



vL,, 



(8) 



In the next sections we will discuss our choices for lim- 
its of integration, normalization wavelength/frequency, 
and SED models. 

5.1. Limits oj Integration 

As was discussed by IMarconi et al.l (J2004D , the full ob- 
served SED includes light that is not directed to the ob- 
server along its original line of sight. Thus, the SED 
determined in this manner is not the intrinsic one. For 
example, most of the IR radiation (~30/um-l/im) is 
believed to be produced in a larg e toroidal dusty re- 
gion b eyond the accretion disk (e .g. jKrolik fc Begelmanl 
119881 lElitzur fc Shlosmanl [2009 ). The radiation from 
this "dusty torus" arises from the reprocessing of higher 
energy photons emitted by the accretion disk and re- 
radiated by the torus in the IR. Any optical-UV radia- 
tion that does not come directly to the observer, but is 
instead reprocessed by dust, is effectively double-counted 
when determining th e true bolomet ri c lurn inosity. Be- 
cause of this effect, IMarconi et al.l (|2004D do not in- 
clude the IR bump when computing their bolometric 
corrections. Their choice of integration limits of 1 /im- 
500 keV is well justified, although strictly speaking the 
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(c) C IV dependent SEDs 

Fig. 14. — Mean SEDs for quasars in two zones of the C IV blueshift vs. C IV EW plane. Panel (a) shows Zone 1: BS < 600 km 
s-i and EW > 50 A (disk- dominated objects in [Richards eFall [20T!^ . Panel (b) shows Zone 2: BS > 1200 km s'^and EW < 32 A 
(wind-dominated objects in [Richards et al.1120111 '1. Panel (c) shows the two curves overplotted. See Figure |8] for explanation of scatter 
points in panels (a and b) and the bottom histogram of panel (c). The z one 1 mean SED shows a larger Balmer continuum (A < 3000A) 
and stronger Lya emission, consistent with eigenvector 1 analyses. Figure [131 shows an expanded version of panel (c). 



optical/UV bump must also be corrected for dust red- 
dening along the direct line of sight; otherwise, the ob- 
served luminosity will be smaller than the true intrinsic 
luminosity. 

Here we describe a similar effect in the hard X-ray band 
(energies larger than 2keV). While it is thought that IR 
flux can include reprocessed disk emission, so too can the 
hard X-ray emission. The accretion disk itself can emit 
thermal soft X-ray photons from the inner part of the ac- 
cretion disk, but hard X-ray photons are believed to come 
from Compton upscattering of accretion disk photons off 
hot electrons in the so-called "corona." Depending on the 
geometry of the corona, it may or may not be appropriate 
to include the hard X-ray part of the SE D in the bolo- 
metric luminosity in the same way that iMarconi et alj 



(|2004D suggest avoiding the IR emission. If the corona 
is a hot spherica l regio n surrounding the black hole 
(jSobolewska et al.l l2004aD then one should exclude the 
highest energy photons from the tabulation of the "in- 
trinsic" luminosity. However, if the cor ona is more like a 
patchy skin to the accretion disk as in ISobolewska et all 
^004b), including the hard X-ray flux in the SED would 
be more appropriate. Our solution here is an agnostic 
one; instead of presenting the full integrated luminosity, 
we will give it in four pieces: 30 /xm-l /im, 1 fiin-2 keV, 
2keV-10keV, and 10keV-500keV (see Table [3]), allow- 
ing the user to determine which approach is best. Al- 
though, given the small amount of data and the large 
uncertainties in the X-ray, we do not recommend tak- 
ing BCs based on monochromatic X-ray luminosities as 
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proposed by iVasudevan fc Fabianl ()2007l ) . 

5.2. Normalization Wavelength 

Typically, bolometric corrections are computed with 
respect to 5100 A and we will report those values herein 
for backward compatibility with previous work. How- 
ever, we will generally report BCs relative to 2500 A. 
There are a number of reasons for this choice. While 
5100 A makes sense for low-redshift sources, especially 
when using the H/3 line to estimate black hole masses 
and then the Eddington ratio from the ratio of Lboi to 
Mqh, that rest- frame wavelength is inaccessible for the 
vast majority of SDSS quasars. The SDSS quasar sample 
peaks at z ~ 1.5, which corresponds to rest- frame spec- 
tral coverage of ^ 1500-3700 A. As such, the majority of 
SDSS quasars have observed flux density measurements 
at 2500 A. For this reason, Richards ct al. (2006b) chose 
to iiT-correct to z = 2 where the SDSS i-band roughly 
corresponds to AoS — 2500 A. Moreover, 2500 A is typ- 
ically used as the optical anchor point in the L^v-ciok 
relationship. In addition to this, the 5100 A luminosity 
always has a larger host galaxy contamination. As such, 
we have chosen to use 2500 A as our fiducial wavelength. 

5.3. Integrated Luminosities and Bolometric 
Corrections 

For comparison with previous work, we compute the 
bolometric luminosity in a number of ways, including the 
tabulation of integrated optical and IR luminosities. In 
terms of bolometric corrections, we start by computing 
BCs for each individual object in our sample, and, with 
those, compute mean BCs for the full sample. In comput- 
ing BCs we used SEDs constructed as discussed above: 
namely, using broadband observations where available, 
the Luv~CKox relationship, and "gap-filling" (with the 
mid-luminosity mean SED) as needed. The overall mean 

BC ' and corresponding standard deviation values 
2500 A "^ " 

are 2.75 ± 0.40 using limits of 1 /im and 2keV (here- 
after called idisk, which avoids the issue of IR and hard 
X-ray double counting); at 5100 A this correspon ds to 
4.33 zfc 1.29. Fo r compa rison with lElvis et all ()1994[ ) and 
[Richards et al.l ()2006al ) we also determine the bolometric 
correction to 5100 A in the range of 30 /im up to lOkcV. 

We find BC ° = 7.79 ± 1.69, as compared to the re- 

.^1 on A ^^^^^ ^^^^^^ 

suits of lElviset all (|1994D and iRichards et"all (|2006a| ) 
who found ll.Stg^a^ ^^^ l^-^ ±2.1 respectively. Using 
the range 1 /^m to 500 keV from lMarconi et all ()2004[ ) we 
find that BC » = 2.97 ± 0.43. 

Recent stud i es by iNemmen fc Brothertoiil ()2010D and 
iRunnoe et al.l (|2012f ) suggest the relationship between 
a quasar's monochromatic and bolometric luminosity is 
nonlinear. As is expected given the nonlinear relation- 
ship used to connect the UV and the X-ray, we also find 
this to be the case; our best fit is a power-law of the 
form: 

log(LBoi) = (0.9869 ±0.0003) log (.^^^^^L^^^^^) 

+ (1.051 ±0.014) (9) 

Note that while we allow for a nonlinear slope, the best- 
fit slope is very close to linear. The errors quoted on this 




Fig. 15. 



Disk bolometric correction as a function of UV lu- 



minosity, L ' . The blue contours and scatter points shows a 

density plot for our quasar sample. The contours indicate the lin- 
ear density of scatter points on the plot. The black dotted li ne is 
the bolometric correction for the overall mean SED (see Figure [TOt 
without taking into account the nonlinear correlation between Luv 
and Lx- The blue dashed line shows the L-dependcnce of the 
BCs assuming the Luv — L^ relation (Steffcn ct al. 2006); this line 
tracks the data points since the individual BCs are also dependent 
on -Luv • The black dashed line shows the same relation with an 
extra bump added in the EUV (see Section |6ll . For comparison 
the relationship from Marc oiu et al.l ([2004) is shown as the dotted 
maroon line scaled to the limits of integration and to 2500A. 



fit are small because of our large sample size and only 
statistical, not systematic, errors are included. 

In addition to mean BCs, we give BCs and error es- 
timates for each of the objects in our sample; these are 
shown by the points and contours in Figure ITS] As did 
iMarconi et al.l (|2004l ). we investigate how the BCs are de- 
pendent on quasar luminosity and, in turn, on the Luv- 
aox relation that was assumed. The black dotted line in 
Figure [15] shows the luminosity-independent BC that one 
gets for our mean SED and assuming a linear dependence 
between Luv and Lx(i.e. re- normalizing the mean SED 
without changing aox). To illustrate the importance of 
the observed nonlinear correlation between Luv and Lx 
we do the following. We re-normalize our mean SED to 
L = values that span the range shown in Figure [TSl 
stepping in small increments. We then truncate the mean 
SED at 1216 A and connect it to 0.2 keV using the Luv- 
aox relationship and assuming a^ — —I beyond 0.2 keV. 
We then calculate the resulting BCs over a range of lumi- 
nosities; the results are given by the dashed blue line in 
Figure [m This line tracks the outliers well, but the core 
of our sample shows a weaker luminosity dependence to 
the BCs (as expected from Equation ^). This weaker 
luminosity dependence may reflect the trend of redder 
optical continua in low-luminosity quasars are bluer op- 
tical continua in high-luminosity quasars, which would 
counter-act the trends in aox with Luv 



We have also included the line found by Marconi et al. 
(|2004| maroon) for comparison. The .Marconi et al 
(2004) line has been adjusted to match our limits of 
integration and the BC normalization wavelength. We 
see that the form of the Marconi line and ours are sim- 
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TABLE 3 
BoLOMETRic Luminosities and Bolometric Corrections 



SDSS ID 



^°g(^^2500A) 



log (h'L » ) 

^'- 5100 A-^ 



BCh 



log(i30Mm-lfim) log (-Lif,ni-2kcv) log (i2kcV-10kcv) log (ilOkcV-SOOkcv) 



587745539970630036 


46.06 


0.03 


45.77 


0.07 


3.08 


0.25 


46.27 


0.03 


46.54 


0.01 


44.75 


0.17 


45.11 


0.11 


587732772647993547 


45.35 


0.07 


44.98 


0.04 


2.59 


0.43 


45.52 


0.02 


45.76 


0.02 


44.24 


0.17 


44.6 


0.11 


587735666930941958 


45.68 


0.08 


45.37 


0.06 


2.83 


0.51 


46.13 


0.02 


46.13 


0.02 


44.48 


0.17 


44.84 


0.11 


587724198277808231 


46.31 


0.06 


46.39 


0.08 


2.31 


0.33 


46.78 


0.02 


46.67 


0.01 


44.97 


0.17 


45.33 


0.11 


587729159520452814 


45.85 


0.05 


45.62 


0.06 


3.14 


0.36 


46.22 


0.02 


46.35 


0.02 


44.62 


0.17 


44.98 


0.11 


588848900997054741 


44.71 


0.06 


44.53 


0.07 


3.62 


0.51 


45.02 


0.02 


45.27 


0.02 


43.79 


0.17 


44.15 


0.11 


588017978340999239 


46.91 


0.06 


46.63 


0.03 


2.49 


0.33 


47.25 


0.02 


47.31 


0.01 


45.38 


0.18 


45.74 


0.12 


587722984428798040 


46.29 


0.05 


46.07 


0.07 


3.33 


0.37 


46.47 


0.02 


46.81 


0.01 


44.92 


0.17 


45.28 


0.11 


587732483289055296 


45.5 


0.06 


45.37 


0.04 


2.2 


0.3 


45.94 


0.01 


45.84 


0.01 


44.31 


0.17 


44.67 


0.11 


587733081880789140 


45.54 


0.05 


45.3 


0.04 


2.64 


0.34 


45.91 


0.02 


45.96 


0.02 


44.39 


0.17 


44.76 


0.11 



Note. — All luminosities are reported in log (erg s ^). This table is published in its entirety in the electronic edition of the online journal. A 
portion is shown here for guidance regarding its form and content. 
'^ Taken over the range of 1 iim to 2 keV. 



ilar, and the offset is due to the fact that the Marconi 
SEDs use a shallower Luv~ctox relation that is extrapo- 
lated to 2 keV (which would predict a soft X-ray excess in 
all quasars), whereas we use a steeper Luv-ckox relation 
extrapolated to 0.2 keV. 

The general shape of the distribution of individual BCs 
shows that the main luminosity dependence of our BCs 
is indeed due to the Luv-aox relation. We do, however, 
see a large spread of ~1 on either side of the dashed line. 
In particular, it is possible for a high-L quasar to have a 
EC that is higher than a low-L quasar even though the 
general trend is in the other direction. 

In order to facilitate the determination of bolomet- 
ric corrections using different limits of integration and 
normalizing wavelengths, in Table [3] we have tabu- 
lated the integrated luminosity and errors for individual 
SEDs over 30^m-l^m, l^m-2keV, 2keV-10keV, and 
10keV-500keV. The second range is our recommended 
range as it corresponds to Ldisk, but t he sum of the 
first th ree ranges matches that used by [Richards et al.l 
(|2006af) and the sum of the last three ranges is that 
used bv iMarconi et al.l (|2004[ ). We further give L °, 

L ', and the resulting EC ° (relative to idisk)- 

2500 A' ° 2500 A ^ ^ _ ' 

We have not corrected for non-isotropic emission (i.e. 
the emitted light is not the same in all directions) in 
our tabulation s : how ever, taking anisotropy into account, 
iRunnoe et al.l (|2012j ) suggest scaling the bolometric lu- 
minosities by 0.75 when calculating the bolometric lumi- 
nosity over the range of 1 /im-lOkeV. 

With individual ECs for each quasar, we are able to 
construct EC-dependent mean SEDs using ECdisk- For 
this we split our sample into four equally populated bins: 
EC< 2.65, 2.65 <EC< 2.87, 2.87 <EC< 3.15, and 
EC> 3.15, each containing ~ 29,000 quasars. Unlike 
the luminosity- and C IV-dependent SEDs, these mean 
SEDs do not truncate the individual SEDs and connect 
the mean to the X-ray, but instead connects each individ- 
ual SED to the X-ray and then takes the mean. This is 
done since each SED needs to be connected to the X-ray 
in order to calculate the EC. Figure [TBI shows the result- 
ing means and Figure 113r d) shows the UV region when 
the SEDs are normalized at 2500 A. From these figures 
we see that the mean SEDs only differ significantly in the 
UV (A < 2000 A) and the largely unknown regions of the 
SED in FUV could have significant effects as compared 
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Fig. 16. — Mean SEDs for four bolometric correction bins. The 
higher BC SEDs show more emission in the UV-EUV part of the 
SED (A < 2000 A) and an overall higher luminosity. The steep 
drop-off of the low-BC SED shortward Lya coul d be caused by 
unacc ount ed for dust reddening (see end of Section [33] for details). 
Figure [T3l shows an expanded version of the UV section of this plot. 

to the differences seen over 1000-2000A where the SED 
is well-measured. Although they were not included in 
the mean SEDs, the 11,468 quasars that show signs of 
significant dust reddening {A{g — i) > 0.3) all fall in the 
two lowest EC bins; this is expected since quasars with 
heavy dust reddening will appear to have a smaller EC. 
This cut only accounts for strong dust reddening; if there 
is only a small amount of dust reddening then a high-EC 
quasar could easily fall into the low-EC bin. We might 
expect the low-EC bin to be contaminated by quasars 
with mild dust reddening, which could be the cause of 
the steep drop-off in the low-EC SED just past Lya. 

6. DISCUSSION: THE UNSEEN EUV CONTINUUM 

One of the reasons that we explored the mean SED in 
subsamples as a function of various quasar parameters 
is that we expect that the diversity of quasar proper- 
ties (especially in broad emission lines) is a direct con- 
seque nce of the diversity in SEDs (e.g., [Richards et al.l 
[201 If ). In particular, quasars with significant UV lu- 
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minosity relative to the X-ray are capable of driving 
strong winds th rough radiation pressure on resonanc e 
line transitions (jMurrav et all 119951: iProgaet al1l2000[ ). 
Give n the nonlinear rel ationship between L^v and Lx 
(e.g., Steffen et al. 2006), one might expect to see sig- 
nificant differ ences in the SEDs of q uasars as a function 
of luminosity. iRichards et al.l ()2011[ ) have argued that a 
potential wind diagnostic may be the properties of the 
C IV e mission line in the context of the lBoroson fc GreenI 
(|1992| ) eigenvector 1 parameter space which looked at 
the difference between quasars using principal compo- 
nent analysis. This was the main motivation for explor- 
ing the mean SED as a function of these parameters in 
Section ID While we saw significant differences in the 
mean SEDs as a function of luminosity, we saw fewer dif- 
ferences between the C IV composites than might have 
been expected. 

While it is possible that the SEDs of quasar ex- 
trema (in UV emission) truly are similar, it is difficult 
to understand how BELR properties could be so dif- 
ferent, yet have similar SEDs. For example, there is 
a large range in the strength of the He II 1640 line 
(with an i onization poten t ial of 54.4 eV) in Figures 11 
and 12 of IRichards et al.l (120111). As such, we are le d 
to a similar conclu si on as [N etzer fc Davidsoiil (1 19 791) . 
iKorista et all (|1997[ ). iDone e t al. (2012I). and iLawrencd 
(|2012f) — namely that the EUV SEDs may be quite dif- 
ferent from the standard power-law parameterization be- 
tween the optical and X-ray. This variation could be in 
a number of forms, including a second bump in the EUV 
or simply that the SED seen by the BELR is different 
from what we view f e.g. . IKorista et aLlll997[ ). 

Here we suggest that the solution may be more than 
just a different mean SED in the EUV than is normally 
assumed, but rather that the shape of the EUV SED 
must be quite different for quasars at opposite extrema 
in terms of their BELR properties; i.e., there is no uni- 
versal quasar SED (even after accou nting for the L^v-ctox 
relationship)! While ' Grupe et al.l ()2010h find that only 
a factor of a few difference can be hidden in the EUV, 
their sample is restricted to low-redshift AGNs, whereas 
our suggestion is that the differences might be signifi- 
cant when considering the extremes of the distribution 
as spanned by the full SDSS quasar sample. In terms 
of a model where the BELR emission comes from both 
disk and w ind components (jCollin-Souffrin et al.l 119881 : 
iLeighlv eTa l. 2004), quasars with strong winds would 
have a weaker EUV conti nuum than quasars wi th strong 
disk components. What IRichards et al.l (|201lD refer to 
as disk-dominated objects have emission line features 
that are consistent with a much har der EUV SE D than 
the wind-dominated objects (Kruczc k et al.ll2011h . If the 
SEDs over the observable range of the EM spectrum are 
similar, the unseen part of the EM spectrum may yield 
very different SEDs. This hypothesis is consistent with 
radiation line driving being sens itive to the UV to ioniz- 
ing flux ratio ()Proga et al .ll20000. 

If this hypothesis is correct, it would have important 
consequences for the determination of bolometric correc- 
tions (and thus L/L^dd ex Lboi/Mbk)- A full investi- 
gation into the range of EUV continuum properties is 
beyond the scope of this paper; however, herein we have 
created some examples showing how different these val- 
ues might be for different assumptions of the EUV SED 
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Fig. 17. — The three models we used to represent the EUV part 
of our SEDs at three luminosities. The Casebeer et al. model is 
shown in red and the Scott+Luv-Lx model in green. The black 
dashed line shows a blackbody with a peak energy of 23 eV and 
normalized so it peaks with log{uL,j) = 44.5. These models show 
the differences between the observed and theoretical models short- 
ward of Ly« as a function of luminosity and illustrate the relative 
level of our hypothetical extra EUV component. 

at the extrema of BELR properties. 

The most obvious deviation from our baseline pre- 
scription (using Luv-ttox to describe the unseen part o f 
the SED) follows from the work of iScott et aD ()2004[ ) 
in the FUV. The UV spectral inde x just shortward o f 
Lyg was found to be different by iTelfer et al.l ()2002l ) 
and IScott et al.l ()2004[ ) (apuv = —1.42 as compared to 
aox= —1.54 for the average quasar luminosity). These 
differences can be explained by a lu minosity-dependent 
spectral index given in Equation ([T]) (jScott et al.ll2004[ ) 
that is similar to that seen for Oox and by the fact that 
these two samples probe different ranges of L^v Thus, 
we create L-dependent model SEDs that do not just con- 
nect 2500 A directly to 0.1 keV, but instead that connect 
2500 A first to 500 A following Equation [1] and then con- 
nect 500 A to 0.2 keV using the L^v-aox relationship to 
set the X-ray continuum level. Such an SED is illustrated 
(green lines) in Figure [T7] for three different luminosi- 
ties, where it can be seen that these SEDs have a small 
(luminosity-dependent) excess of EUV flux as compared 
to a power-law fit between the optical and the X-ray. 

We compare this to model SEDs taken from 
iCasebeer et al.l (I2006D that are adapted from the 
CLOUDY (lF^rland||2003) so-called "AGN" continuum. 
The functional form of this continuum is given in 
[Casebeer et al., (|2006l Equation (Al)). Essentially this 
SED consists of a power law representing the optical/UV 
continuum with exponential cutoffs in the infrared and 
UV, plus a power law in the X-rays, where the normal- 
izatio n of the two c o mponents is set by ctox (calibrated to 
Wilke s et allll994[) . ICasebeer et"all ((2006D used a range 
of Q^ox and UV cutoffs to test the influence of the SED 
on emission line ratios (and not to determine BCs); we 
show in Figure [17] the ones with kTcut = 50, 37, 27 eV 
and Qfox = —1.5, —1.57, —1.66, from lowest to highest lu- 
minosity respectively (chosen to have the same optical 
luminosities as the Scott lines). 

While we cannot directly measure the EUV part of the 



17 



1000 A 



100 A 



10 A 




7000 



6000 



I 5000 



BC=2.75±0.80 
BC=2.86±0.91 
BC=5.05±1.95 
BC=3.12±1.01 



O 



4000 



ou 3000 



■ 2000 



1000 



"15.0 



15.5 



16.0 16.5 

\og(v[Hz]) 

(a) EUV models 



17.0 




2 4 6 

Bolometric Correction at A=2500A 

(b) Bolometric Corrections 



10 



Fig. 18. — (a) Mean SEDs for four different models: Luv— iJx from lSteffen et aTl ||2006|. blue), the Luv—L^ relation with a 23 eV blackbody 
with a peak value of log(i'L^) = 45 (black), using the EUV extension from |Scott et al.l 12005, green), and the Casebeer et al. model 



-Lx from lSteffen et al.| ll200g. blue), the Lu 

ilUV extension from jScott et al.l pOQl, gr 
dUascbccr ct al. 2006, red), (b) Bolometric corrections for the four models in (a) where the bolometric luminosity is taken from 1 ^m to 
2keV (idisk)- The errors in the legend give the 2(t levels of the BC distributions. Models similar to all of those shown here are arguably 
well-justified and used in the literature. Their differences emphasize the need to better understand the distribution of EUV continua in 
quasars. 



spectrum, it is interesting to compare the Scott and Case- 
beer et al. SEDs in the '^800 A range. For low-luminosity 
quasars, the data-driven Scott-based SED and the more 
theoretical (but empirically calibrated) Casebeer et al. 
SED (chosen to have the same ttox), follow the same up- 
ward trend. This leads one to wonder if low-luminosity 
objects may indeed have much more EUV flux short- 
ward of 500 A (i.e., that the Casebeer et al. 2006 model 
is correct) than is typically assumed when applying the 
Luv~aox prescription, as the EUV is unconstrained for 
the green curves between 500 A and 0.2 keV. However, 
for the high-luminosity quasars, the Scott and Case- 
beer et al. SEDs do not follow the same trend e ven at 
~800 A, which suggests a model like the Casebee r et al.l 
([2006) SED likely overestimates the EUV emission in 
the most luminous quasars. Indeed, a comparison of the 
UV emission lines in the extrema shown in Figure 12 of 
[Richards et al.l (|2011D suggests that the SEDs of quasars 
with "disk"- and "wind" -dominated BELRs may be very 
different in the range of ^50 eV whe re the ionization po- 
tentials of many UV lines lie (e.g., Richards et al.|[201lL 
Figure 13). 

Since the wind-dominated objects are more luminous 
on average than the disk-dominated objects (as can be 
seen in Figure [HT c)). we also consider a model where 
we add a blackbody component of fixed luminosity, such 
that it is significant in low-luminosity sources, but con- 
tributes little to the EUV continuum for high-luminosity 
sources. 

To construct a toy model to illustrate such a situa- 
tion, w e start by ext e nding our SEDs to the X-ray us- 
ing the 'Steffe n et al.l ([2006) L^^-Uox relation and then 
added a blackbody peaking at 23 eV. This peak was cho- 
sen in order to add photons of the typical energy needed 
to produce the spectral lines we see (e.g., C IV and 
He II). The normalization of the blackbody was taken 
to have log {vL^) = 44.5 at its peak (see the dashed 



black curve in Figure [T7|) so that the highest luminos- 
ity quasars remain unaffected but the lower luminosity 
quasars gain an extra feature which has the correct sign 
to explain the broad range of EWs for, e.g., C IV and 
He II ([Richards et al.ll2011[ ). This model leads to a larger 
spread in their bolometric corrections than when no ex- 
tra component is used (see Figure ITSfb)). While we lack 
direct observational evidence for such a component, the 
range of emission line strengths for lines with ionization 
potentials of ^50 eV, suggestion that some, but not all, 
quasars have a need for more ionizing fiux in this region 
of the EM spectrum in order to understand their emis- 
sion line properties. We include it here as a way to illus- 
trate how significantly such a component would change 
the bolometric corrections of low-luminosity quasars. 

We compare these four models, which were chosen to 
bracket the range of reasonable EUV shapes for a given 
ttox, shown in Figure rTSf a) for a single luminosity. Fig- 
ure ITST b) then shows how the BC distribution would 
change under the assumption of each of these models. 
The biggest deviation from the standard model results 
when using the Casebeer et al. SEDs. These SEDs may 
have too much EUV fiux for high-luminosity sources, 
while being more reasonable at low-luminosities (see Fig- 
ure [TT]) . On the other hand, the Scott-based and ex- 
tra EUV component SEDs depict a more subtle shift 
in the mean BCs. More importantly, however, is that 
these deviations from our standard SED produce BC 
changes that are systematic with L^v Although there 
is a lack of data in the EUV, in Figure [151 the dashed 
black line shows how much we would expect the BCs to 
change for low-luminosity sources if there were an extra 
EUV continuum component. Adding the L-dependence 
from Scott to this model would make this contrast even 
stronger. 

While this extra EUV component is speculative, the 
ramifications are such that it is important to consider the 
possibility. Moreover, the L-dependence of the FUV con- 
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tinuum based on the work of iScott et al.l ()2004D already 
suggests that a correction in this direction is needed. 
Achieving a better understanding of the EUV contin- 
uum is clearly of importance for understanding quasars 
physics, as a change of a factor of a few in BCs for low- 
luminosity sources translates to the same correction fac- 
tor for both Lboi and the Eddington ratio (i.e., accretion 
rates). Ideally, we would also like to consider the BC 
distribution as a function of quantities like the accretion 
rate, but such analysis is difficult without first under- 
standing any systematic effects in the unseen EUV part 
of the spectrum. 

7. CONCLUSIONS 

We have compiled a sample of 119,652 quasars detected 
in the SDSS. These data are supplemented with multi- 
wavelength data spanning from the mid-IR through the 
UV (Table [1]). This data set was used to construct a new 
mean SED consisting of 108,184 non-reddened quasars, 
with rest frame coverage from '^ 20 /zm - 912 A (Ta- 
ble [2|). By splitting our sample into luminosity bins we 
constructed three luminosity-dependent mean SEDs and 
found apuv to be dependent on luminosity in the sense 
that more luminous quasars hav e redder aFuy con tinua 
(Figure [iSlJb)), consistent with IScott et al.l ()2004[ ). In 
addition, the high-luminosity quasars also show signs of 
having bluer optical continua (Figure I13f b)) and more 
hot dust emission than the low-luminosity quasars (Fig- 
ure [THKa)). When splitting our sample based on C IV 
properties we saw differences in the Balmer continuum, 
Lya and C IV (by construction) line strengths (Fig- 
ure [T3l^c)) that are consistent with eigenvector 1 trends 



(jBrotherton fc Francislll999f) . 

We also constructed SEDs for each quasar and, from 
those, found bolometric corrections (Table [S]). The over- 
all mean is BC 



2500 A 



2.75 ± 0.40 using integration lim- 
its of 1 fim- 2 keV. While the range of bolometric cor- 
rections indicated by the IR through NUV data alone 
is reasonably small, there can be significant changes in 
the distribution of bolometric corrections when different 
models are assumed in the unseen EUV part of the SED 
(Figure [TSTbll. Although more work has to be done to 
determine which model should be used in the EUV, it 
is nevertheless clear that it is important to consider po- 
tentially significant differences in the EUV part of the 
SED at the extrema of quasar properties (e.g., in lumi- 
nosity and in eigenvector 1 parameter space). In future 
work, we will further explore the cause(s) of the width in 
the bolometric correction distribution (after correcting 
for known luminosity trends) and will characterize the 
systematic error in the accretion rate estimates. 
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